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'"^ Abstract. We use sub-parsec resolution hydrodynamic resimulations of a Milky 

' Way (MW) like galaxy at high redshift to investigate the formation of the MW 

satellite galaxies. More specifically, we assess the impact of supernova feed- 
back on the dwarf progenitors of these satellite, and the efficiency of a simple 
instantaneous reionisation scenario in suppressing star formation at the low- 
\^ mass end of this dwarf distribution. Identifying galaxies in our high redshift 

(-H simulation and tracking them to z=0 using a dark matter halo merger tree, we 

^ I compare our results to present-day observations and determine the epoch at 

which we deem satellite galaxy formation must be completed. We find that only 
the low-mass end of the population of luminous subhalos of the Milky-Way like 
galaxy is not complete before redshift 8, and that although supernovae feedback 
^ reduces the stellar mass of the low-mass subhalos (M< IO'Mq ), the number 

—I of surviving satellites around the Milky- Way like galaxy at z = is the same 

in the run with or without supernova feedback. If a luminous halo is able to 
avoid accretion by the Milky- Way progenitor before redshift 3, then it is likely 
to survive as a MW satellite to redshift 0. 
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1 Introduction 

The "Missing Satellite Problem" arose from the discovery that pure N-body simulations 
of MW-sized halos contain at least an order of magnitude more dark substructures than 
observed luminous satellites (fT7llT3l l. Although this discrepancy has recently been allevi- 
ated by the detection of a host of of new MW satellites (14, 29 ), a non negliglible (factor 
of a few) overabundance of dark substructures remains. Processes such as reionisation 
(|6l [TT] m |251 ) and supernova feedback (|2l [3j have been invoked to reduce the number 
of satellite galaxies by heating and expelling potentially star-forming gas from low-mass 
halos. 

To test this claim, we are running extremely high resolution hydrodynamic cosmolog- 
ical resimulations of a MW-like galaxy and its surroundings, improving on the work of 
other authors ([T9ll3Tb by at least an order of magnitude in physical resolution. 
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Fig. 1. Comparison of stellar masses of galaxies at z=8.5 in the Cooling and Feedback run whose 
particles end up inside the MW virial radius at z=0, showing the fractional difference between 
the values in the two runs vs total virial mass in Adiabatic Run. A positive value indicates that 
the inclusion of supernova feedback enhances star formation, while a negative value means that it 
suppresses star formation. The circle size indicates stellar mass (see legend). Red circles indicate 
halos that survive as MW satellites to z=0; blue circles are completely disrupted by merger events 
by z=0. 



In this paper, we address the issue of whether MW satellite formation was complete by 
the epoch of reionisation or not. In addition, since we run identical simulations, but with 
and without supernova feedback, we assess the role of supernovae during the formation 
of satellite galaxies at high redshift. Since our resimulations containing star formation 
do not reach z=0, we compare our results to observations of MW satellites by twinning 
the halos containing galaxies in these high-redshift simulations with halos in an adiabatic 
gas simulation run to z-0. Along the lines of 16_, 26^ we track dark matter halos to z=0 
to determine the properties of surviving halos. We then proceed to quantify the effect of 
various physical processes on dwarf galaxy formation at high redshift, and determine the 
importance of such processes in shaping the main properties of the MW satellite galaxies 
we see today. 



2 Numerical Simulations 

In this work, we analyse three resimulations of a Milky Way-like galaxy from the SS=ij!sll 
(Nut) suite of simulations (Powell:2010inprep, in prep) named "Adiabatic", "Cooling" 
and "Feedback". We carry out these simulations using the Adaptive Mesh Refinement 
code RAMSES (l28T l. All runs contain dark matter and baryonic gas. For the Cooling 
and Feedback runs, we include metal dependent radiative energy losses by the gas and its 
heating by a spatially-uniform UV background (12) after z=8.5, along with star forma- 
tion. The threshold density for star formation is 10^ atoms/cm^ with an efficiency of 1%. 
For the Feedback run, we include supernova feedback, which also enriches the gas with 
metals. These simulations have a maximum physical resolution of 0.5pc, with a minimum 
mass of 5.6xlO'*M0 for dark matter particles and 76Mo for star particles JTOb . With this 
mass resolution we are able to capture individual supernovae, although we still treat star 
particles as populations of stars rather than individual stars. 
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Fig. 2. Star formation histories in the Cooling run (left) and Feedback run (right). The vertical 
line at z=8.5 (look-back time 13.065Gyr) shows the point at which reionisation is turned on in the 
simulations. The grey region shows the detectable star formation rates as determined bvl32l 

3 Tracking Halos 

We identify dark matter (sub)halos and clusters of stars using HaloMaker with the Most 
Massive Subhalo Method ( [30] l. We then associate each cluster of stars to its nearest dark 
matter (sub)halo. For the adiabatic run, we use TreeMaker (l30b to build a halo merger tree 
of the MW-like halo from our initial conditions to z-Q. We then attach all dark matter 
(sub)halos in the Cooling and Feedback runs to their nearest "twins" in the Adiabatic 
run (based on the number of dark matter particles shared between the two halos) for a 
selection of redshifts before and after reionisation between z=ll and z=8 ( 10). 

Using the Adiabatic run merger tree, we track each "twinned" halo to z=0 and find 
whether it has survived as an independent halo, as a subhalo of the MW, or has been 
completely disrupted in a merger with the MW or another halo. We posit that since all 
our galaxies are embedded within a dark matter (sub)halo, if this host halo survives to 
z-Q then so does the embedded galaxy. 

From the 87 galaxies formed in the Feedback run at z=8.5 (just prior to reionisation) 
we find 13 galaxies survive as MW satellites to z=0, with 4 being destroyed in mergers 
with other satellite progenitor galaxies (6 in the Cooling run) and 7 1 destroyed in mergers 
with the MW (69 in the Cooling run). At z=8 (51 Myr after reionisation), one further 
satellite galaxy progenitor is formed that survives to z=0. 

4 Feedback and Reionisation 

Comparing the Cooling and Feedback runs, we find that galaxies with total mass under 
IO^Mq generally exhibit suppressed star formation when our model of supernova feed- 
back is included. By contrast, galaxies with total mass over lO^M© have similar masses in 
both runs, with some galaxies showing signs of enhanced star formation when supernova 
feedback is included (figure[T]l. 

The latter is a consequence of enhanced cooling and hence higher star formation rate 
due to the introduction of metals from supernovae into the interstellar medium. However, 
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Fig. 3. Cumulative velocity dispersion distributions of halos with galaxies and/or star clusters em- 
bedded in the Cooling (Blue) and Feedback (Red) runs, tracked to z=0 as explained in section |3] 
Observations (black line) use data from ( I15| [ni2ll33t. The simulated data is given as the maximum 
circular velocity / VS. The error bars represent the sampling error on the data in each bin. 

supernova feedback does not appear to affect the number of galaxies formed, with roughly 
similar numbers of galaxies formed in both the Cooling and Feedback run. Further, the 
same MW satellites are found in both the Cooling and Feedback run when tracked to z=0 
as described in section [3] 

We find that reionisation does not immediately halt star formation, even in the small- 
est halos, in agreement with '3T (see figure [2|i. There appears to be a slight downward 
trend in the star formation rate after reionisation, however. As the simulations continue 
beyond 200Myr after reionisation, we will monitor whether this trend continues. Since 
we implement a uniform ionising background, our galaxies cannot self-shield. We plan 
to address this issue in future studies by including radiative transfer of hydrogen-ionising 
radiation. 

5 MW Satellites Today 

As described in section |3] we track the halos containing galaxies in the Cooling and 
Feedback runs to z=0 using the Adiabatic run merger tree. We find that the most important 
criterion for a halo to survive to z=0 is its redshift of capture, with halos captured before 
z=3 being completely disrupted by z=0. However, this redshift value may be dependent 
on the host mass. 

In figure |3] we compare the observed stellar velocity dispersions to our simulated 
galaxies' dark matter velocity dispersions at z=0 in the Adiabatic run, since estimates 
of virial mass from observations have typically large errors OJ' For all sampled redshifts 
(z=l 1 to z=8), we find too few satellites when compared with the observed MW satellites, 
with better agreement if we wait until lower redshift before sampling the galaxies. As- 
suming this trend continues, this suggests that satellite galaxy formation is not complete 
by reionisation. Indeed, we find one galaxy formed between z=8.5 and z=8, after our 
implemented epoch of reionisation. Furthermore. fl4l l29l suggest that the current sample 
of observed MW satellites is incomplete. Since our runs containing star formation have 
not reached z=0, we cannot comment on the completeness of our simulations. 
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6 Conclusions 

We find that the velocity dispersions of satellites we produce underproduces those of ob- 
served Milky Way satellites, suggesting that satellite progenitor formation is not complete 
by reionisation. Star and galaxy formation continue past reionisation, with little sign of 
reionisation immediately halting the star formation rate. Supernova feedback reduces the 
star formation in halos with total mass < lO^M© , but not the number of galaxies formed. 
Finally, we find only halos captured after z=3 survive to z=0 as Milky Way satellites. 
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